Scaling of success probabilities for linear optics gates 
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By using the abstract linear-optical network derived in [S. Scheel and N. Liitkenhaus, New J. 
Phys. 6, 51 (2004)] we show that for the lowest possible ancilla photon numbers the probability of 
success of realizing a (single-shot) generalized nonlinear sign shift gate on an (N + l)-dimensional 
signal state scales as 1/N 2 . We limit ourselves to single-shot gates without conditional feed-forward. 
We derive our results by using determinants of Vandermonde-type over a polynomial basis which is 
closely related to the well-known Jacobi polynomials. 

PACS numbers: 03.67.-a,42.50.-p,42.50.Ct 



I. INTRODUCTION 

One of the promising routes to implementing small- 
scale quantum information processing networks for tasks 
in quantum communication or cryptography is by us- 
ing linear-optical networks. Already a number of exper- 
iments have been performed that accomplish the task of 
generating controlled-NOT gates as the basic nontrivial 
two-mode operation 0, Q • The necessary nonlinearities 
on the single-photon level are produced via conditional 
measurements 0,0. 

Some two-mode gates such as the controlled-o^ gate 
can be generated by acting separately on both modes 
within a Mach-Zehnder interferometric setup. In that 
way the complexity of a two-mode gate is reduced to a 
single-mode gate, the nonlinear sign shift gate. A number 
of recent theoretical works 0, 0,0, give a variety of 
networks capable of performing the nonlinear sign shift 
gate with the following transformation rule: 

c |0) + ci|l) + c 2 |2) ^ Co |0) + Cl \l) - c 2 |2) (1) 

which can serve as a building block for the two-mode 
controlled-phase gate. The crucial question for the abil- 
ity to concatenate linear-optical gates are their proba- 
bilities of success. Recently, it has been shown that the 
nonlinear sign shift gate cannot be realized with a suc- 
cess probability of more than 1/4 if one abstains from 
using conditional feed- forward [9j, llfJl which is stronger 
than the bound of 1/2 obtained in [Tl|. 

These results were obtained by considering an abstract 
type of network that includes all possible networks by 
dividing it into an 'active' beam splitter 'A' that mixes 
signal and ancilla states, a preparation stage 'P' and a de- 
tection stage 'D' that are there to generate and disentan- 
gle a purpose-specific ancilla, respectively. The general 
scheme is reproduced in Fig.^ In short, the abstraction 
consists in the observation that each U(A^) network can 
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FIG. 1: Preparation stage 'P' of the ancilla state that is fed 
into the 'active' beam splitter 'A', and decomposition stage 
'D' for detection. 



be decomposed into a triangle-shaped network of at most 
N(N — 1)J2 beam splitters (and some additional phase 
shifters) [l^] which in fact is nothing else but doing a 
Householder-type diagonalization of the unitary matrix 
associated with the network. Then, it is easy to see that 
the signal mode can be chosen to impinge only on a single 
beam splitter (here called 'A') with the rest of the net- 
work conveniently divided into the ancilla preparation 
and detection stages 'P' and 'D', respectively. 

What we will try to do in this article is to answer the 
question how the probability of success scales as the sig- 
nal state becomes higher-dimensional, thereby giving a 
hint on the possible scaling law for multi-mode quantum 
gates (by virtue of the above-mentioned fact that general- 
ized nonlinear sign shift gates can be combined in Mach- 
Zehnder interferometer-type set-ups to yield multi-mode 
controlled-phase gates). However, as we have seen pre- 
viously, a general answer is hard to obtain. Fortunately, 
it turned out that the most important special case can 
be dealt with analytically. It has been argued already 
111 and confirmed in |l(j that the success probability, 
taken as a function of the chosen ancilla photon num- 
bers, reaches its maximum if the ancilla contains as few 
photons as possible. In addition, the restriction to the 
lowest possible photon number is justified by the fact 
that in this case the least quantum-state engineering is 
needed to generate the ancilla in the preparation stage 
'P' and decoherence affects the ancilla only minimally. 

This article is organized as follows. In Sec. [H] we pro- 
vide the necessary framework that is suitable for our set 
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task. Specializing to the lowest admissible photon num- 
bers, we show in Sec. 11111 that there exist networks that 
realize the generalized nonlinear sign shift (to be defined 
in Sec. [HJ with a probability of success of 1/N 2 which 
constitutes the main result of our paper. Some details of 
the calculation are provided in the Appendices. Conclu- 
sions are drawn in Sec. IIVI 



II. 



(N + l)-DIMENSIONAL SIGNAL STATE 
WITH iV-DIMENSIONAL ANCILLA 



A possible building block for multi-mode quantum 
gates via some Mach-Zehnder interferometer-type set- 
up is the generalized nonlinear sign shift gate. Let us 
therefore consider an (N + l)-dimensional signal state 



JY 



(2) 



k=0 



The generalized nonlinear sign shift gate is defined by the 
following transformation rule: 



cn 



-Cat 



(3) 



leaving all other coefficients untouched. As thoroughly 
described in not all ancilla states are admissible. In 
fact, the ancilla must consist of a superposition of states 
with fixed photon number. Moreover, the number of pho- 
tons that are detected after the disentangling stage 'D' 
in Fig. n must be the same as the initial ancilla photon 
number. The reason for this requirement is essentially 
that the Hilbert spaces of incoming and outgoing signal 
states have to be isomorphic. 

We will restrict ourselves to the case in which we use 
an A-dimensional ancilla for transforming an (N + 1)- 
dimensional signal state. This choice is somewhat restric- 
tive but, in view of minimizing experimental resources, 
the most accessible choice. Let us denote the N different 
photon numbers in the ancilla by ni, I = 1 . . .N. Then, 
the ancilla can be represented in the form 



N 

i=i 



(4) 



with J2i if = 1 (the expansion coefficients can be taken 
to be real) and n — J2i n i- The state \ A n _ ni ) is arbitrary 
as concerns the number of modes as long as it contains 
exactly n — ni photons. The advantage of introducing the 
abstract network in Fig. ^ is apparent since the |A n _ ni ) 
do not play any role in what follows. 

Now we transform the signal and ancilla states at the 
'active' beam splitter 'A' and use the fact that the matrix 
elements (k,ni\U\k,ni) can be written as 



{k,m\U\k,ni) = (T*) 



-fcp^0,n,-fe)^2|j.|2 



(5) 



where P, 



(o,m-fe) 



(x) is a Jacobi polynomial defined by 



P 



(0,nj-Jfe) 



(2\T\ 



E 



\T\ 



iy 



.(6) 

and T the (in general complex) transmission coefficient 
of the beam splitter 'A'. Then, the (unnormalized) states 
\ipk) {k = 1 ... A + 1) after the transformation can be 
defined as 



N 

E 

i=i 



7; (T * )ni -k+l p(0,n, (p) \ ni) \ An _ ni) 



(7) 

where we have abbreviated the argument P — 2\T\ 2 — 1 
of the Jacobi polynomial to emphasize its dependence on 
the permanent of the beam splitter matrix. This is a 
reminder of the fact that matrix elements of unitary op- 
erators can always be written as permanents of a matrix 
associated with the symmetric tensor power of the beam 
splitter matrix [l3| . 

The state we project upon is denoted by 



N 

E< 



\ni)\A n . 



(8) 



where |a/| 2 = 1. Then, in order to realize the general- 
ized nonlinear sign shift, we need to fulfil the N equations 



(VIVw+i) = -(V#*>, k = l...N. 



(9) 



That specifies all weights ai as well as the transmission 
coefficient T which can be used to compute the probabil- 
ity as 



P = |(?/#fc)| 2 for any k. 



(10) 



In fact, Eq. © is a homogeneous matrix equation of the 
form 



(11) 



i=i 

with matrix elements 

(!) , 

o-ki = air + a 



hi 

(T*) n >- 
+(T*) n 



(2) 

:— fc+lp(0,nj-A+l) 



fc-1 



(12) 



From linear algebra it is known that the homogeneous 
system of equations <|llll has nontrivial solutions only if 
the determinant of the coefficient matrix vanishes. That 
gives a condition for the transmission coefficient T . 

From the expansion theorem of determinants we imme- 
diately have the solutions for a/7/ given in terms of the 
cofactors A^i of the matrix elements aki, i-e. ai"/i oc Aki- 
Because the success probability (|1U|) is nothing else but 
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P = I Ei a '7i a fcf| 2 : ^ is clearly maximized if the magni- 
tudes of ai and 7/ are the same. From the normalization 
condition for \tp) we then get the solution 



to be computed. The first term is just det(oL ), and the 
other N terms are obtained by replacing the mth column 
of (cijy ) by the corresponding column of (ajy )• Thus, 



for some k . 



(13) 



™ (EM kl \r 

We will shortly see that the choice of k does not matter. 

III. MAXIMAL SUCCESS PROBABILITY 

In most cases, one needs to refer to numerical methods 
to maximize the success probability l|13[l for a given set 
of ancilla photon numbers nj. However, there is one case 
in which an analytical result can be obtained. For this 
purpose, let us consider the photon numbers n; = I — 
1. This is the ancilla with the lowest possible photon 
numbers and as such the experimentally most important 
situation. Then, we find that there exists at least one 
solution of the form 



T = 1 — 2^ N 



(14) 



For N — 2, i.e. the nonlinear sign shift gate with up 
to two photons in the signal state, this is the solution 
T = 1 — \[2 encountered in 0, 0|. For this particular 
choice, the maximal success probability reaches 



Pn 



1 

TV 2 " 



(15) 



This result is obviously true for N = 1 since the corre- 
sponding gate can indeed be done deterministically, and 
for N = 2 that is the result obtained in |j, [Hj • Note the 
surprising fact that the success probability does not scale 
exponentially but rather quadratically with N. Equa- 
tion (|15fl constitutes the main result of this paper. We 
will spend the following subsections deriving it. 



det(a fe i) 



det(4f) + 



N 

E 



detfa 



(2) . J2) 



A7 



4 fem 



b km> 



(17) 



All higher-order replacements vanish due to the fact that 
determinants of matrices with more than one constant 
column vanish identically. 

(2) 

Let us first consider the term det(a kl ). By inspection 



of the structure of the matrix elements 1 
that this determinant can be written as 



i {2) 
'ki 



det(ag>)=det P™-" (2\T\* 1 



(0,1!-*), 



one realizes 



(18) 



Due to the special choice of photon number n; = I — 
1, all prefactors containing the transmission coefficient 
cancel out. The next observation is crucial for all the 
following. The Jacobi polynomials pjfj~ k) (2\T\ 2 - 1), 
being polynomials of order k — 1 in 2\T\ 2 — 1, are also 
polynomials of order k — 1 in the photon numbers n\ — 
l — l. This can be seen as follows. Recall the definition of 
the Jacobi polynomials [Eq. |JB}, we briefly revert to the 
general notation n; for the photon numbers] and rewrite 
the binomial factor ( n ' +m ) as 



m + m 
rn 



^ m ^ m 



(19) 



p=0 



which itself is a polynomial of degree m in rii . Hence, the 
Jacobi polynomial p^°' n ' fc -'(2|T| 2 — 1) is also a polyno- 
mial of degree k in the photon numbers n; . The cr^]_ p are 
the elementary symmetric polynomials of degree to — p of 
to variables which in our case are just the integers 1 ... to. 
We will now define a new class of polynomials (n) by 
setting 



A. Determinant of the coefficient matrix 

In the calculation of the projection vector we 
encountered the matrix consisting of the elements a^i 
[Eq. (|12l) ]. As one can see, the e£ are actually indepen- 
dent of the row index k and therefore constant columns. 
In addition, each column consists of two sub-columns, 
and the determinant can thus be expanded as a sum 
of determinants, each of which contains one of the sub- 
columns 01- Hence, 



det(a fe; ) = ^2 



det(47< 



(16) 



mi,...,mjv— 1 



Together with the property that a determinant vanishes 
if it contains more than one constant column, this means 
that only N+ 1 out of the total 2 N terms eventually have 



S k M = Ck v n \ = nP i J2 

p—0 p—0 m—p 

These polynomials fulfil the identity 



k \ (x 2 - 1) 



p(°>»<- fe )( 2 |T| 2 -l). 



(m) 
m—p " 

(20) 



(21) 



From their definition, Eq. (|20|) . it is easy to show that 
these polynomials fulfil a three-term recursion relation. 
Some properties of the S'-polynomials are collected in Ap- 
pendix^ 

With this definition, we can finally rewrite the deter- 
minant H18fl with the help of the S'-polynomials as 



det(4 2) ) = det S^l>(m) , 



r(|T|), 



(22) 



which is nothing but a Vandermonde determinant over 
the S'-polynomial basis. For later use, let us denote this 
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Vandermondian by Vs,jv({h;}), where the subscripts S 
and N denote the polynomial basis and the dimension, 
respectively. For such a determinant we have the relation 



(see, e.g. UJ) 



JV-1 

Vs } N({ni}) = V N ({m}) c fefe 

fe=0 



(23) 



where V/v({n;} is the usual Vandermonde determinant of 
the photon numbers rii. This implies that 



det(og>) = 



(24) 



Specializing again to ni = l — l, in which case Vjv({nj}) 
fL(Z — 1)!, we are left with 



det(4f) = (|T| 2 -l)f(^) 



(25) 



Next, in order to compute the determinant in Eq. I|17|) 
we also need the terms in which the mth column of {air ) 
has been replaced by (dl-j). Expanding the determinants 



dct(a 



obtain 

JV 



(2) . „(2) 
kl 



a km) w ith respect to the column m, we 



JV 



E **(«S> : 4 2 > - «£) - E 



2) 
nm J 



(26) 



m— 1 



m,n— 1 



(21 (2) 

where the A n m are the cofactors of the matrix (a kl ). 
They are, disregarding some prefactors, nothing but gen- 
eralized Vandermondians over the S'-polynomial basis 
with one gap, denoted by VsN^Jn-i}' w ith the polyno- 
mials of order n — 1 and photon number n m = m — 1 
missing. With that, Eq. 126(1 becomes 



JV 

E 

m— 1 



det(ag> : a (2) 



km 



53(-l) JV - n (T*)"- 



■JV-1 



n=l 
JV 



E(-l) m 4 TD (m-l)^ (m - ] 



S,JV\{n-l} 



(27) 



m— 1 



The last line in Eq. I|27|l is just the definition of 
Ks,jv+i\{>i-i}- Such a quantity can be computed in the 
following way. Analogously to the elementary symmetric 
polynomials a^' defined in Eq. 1|19|) we define the elemen- 
tary symmetric polynomials Si(x; N; N) with respect to 
the S'-polynomial basis by 

I 2 _ i \JV N N 



These polynomials Si(x; N; N) are ith-order symmetric 
polynomials over the N different photon numbers ni, I = 
1 . . . N. Expanding Vs,n+i on one hand with respect to 
a Vandermondian of lower order and on the other hand 
with respect to its last column, viz. 



(x 2 - 1 ) N 

ivT 



Vs,n+i 



JV 



J\( n N+i - rii)Vi 



S,N 



1=1 



JV 



E(-l) JV " m ^ ) ^JV + i)^,JV+i\ {m} (29) 



m=0 



and using the definition l|28[) for the symmetric polyno- 
mials, we obtain that 

Ks,JV+i\{n-i} = s N - n+1 (\T\;N;N)V s , N (30) 

(see also 0). 

The symmetric polynomials Sj(x; N; N) can be com- 
puted by inserting Eq. (|19|l into Eq. I|28|l . When specify- 
ing yet again to the choice of photon numbers ni =1 — 1, 
the solution is simply 



.s„(|71:.Y:.Yl= (*| T 



2k 



(31) 



which can easily be checked by inserting into Eq. I|28|) and 
using the definition i|20[) for the S-polynomials. Note that 
the symmetry with respect to the interchange of any two 
photon numbers n k = k—1 and ni = Z — 1 is hidden inside 
the binomial coefficient in Eq. (|31|l . For a constructive 
proof of the statement l|31(l and generalizations of it, we 
refer the reader to Appendix iBl 

Combining the knowledge about the symmetric poly- 
nomials and the Vandermondians defined earlier, we ar- 
rive at the expression 



JV 



n,m=l 



i=i 

JV 

= J2 { - 1) N-n {TT -N-ly s ^ {n _ i} 



n=l 
JV 



N 



= ^(-l)^-"(T*)"- JV - 1 ( _ )\T\ 2{N - n+1) V s , N 

n=l ' 

= (|T| 2 - l)-^- 1 [1- (1-T) N ] . (32) 

(2) 

Combined with the result for the determinant det(aL ) 
[Eq. (|25|) ] this yields the final result for the determinant 
of the coefficient matrix as 

det(a kl ) = (|T| 2 - l)^^ [2 - (1 - T) N ] . (33) 

A nontrivial condition for this determinant to vanish is 
thus 



(28) 



T =1 — 2 1 ' N 



(34) 
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just as we set out to show. What we have seen is that cal- 
culating the value of the transmission coefficient T of the 
active beam splitter 'A' is still mathematically involved 
despite the reduction of the problem to the abstract net- 
work. 



B. Success probability 

In order to compute the success probability 1)13(1 . we 
need to calculate the cofactors Aki of the generalized 
Vandermondian (aki)- It turns out that it is useful to 
choose k = N in this relation since the cofactors Am 
are essentially generalized Vandermondians over the S- 
polynomial basis which can be computed in the same 
way as shown above. Although this particular choice for 
k is not unique, it represents the most easily tractable 
situation. Clearly, any other fc could have been chosen. 
However, since the solution T = 1 — 2 1 / A ' is unique and 
is a simple root of Eq. the result cannot depend on 
the particular choice of k. Since now Tel, we can drop 
asterisks to simplify notation. 

The numerator in Eq. I|13|) can be calculated in the 
following way. By using simple row manipulations, it is 
shown to be equal to 



JV-l N 



N 



N—l 

E< 

n=l 



,(2) . (2) 



' k! 



Nl 



{<$ : off h 




a (1) ) 

a Nl> ' 


(35) 



where we have made use of the fact that the determinant 
det(cifcj) vanishes for T — 1 — 2 1 / JV . All subsequent for- 
mulas have to be understood with the solution for T in 
mind. 

(12) 

Inserting the definitions for the matrix elements a kl ' , 
it turns out that this is a Vandermondian with one gap. 
More precisely, 



aei(a k[ . a Nl h-> a Nl ) - 

_ T -1 V (N) 

1 * S,N+l\{N~l} ■ 



(36) 



Expanding this determinant with respect to the last row, 
wc get 



^s,n+i\{n-i} ~ E ( ^) N+k+1 ^N\^Ws,N-i ■ (37) 



fc=0 



The Vandermondian over the polynomial basis can eas- 
ily be computed from the corresponding Vandermondian 
over the power basis as (see Appendix 



(38) 



With that and the definition of the S'-polynomials one 
finds that 



N-l 



k=0 



E(-i) mi r fc _ V(*) 



N 



k=0 m=l 
N 

= E 

771 — 1 



N-l 
k 



k + m 
m 



N 



(T 2 - 1)' 



(T - l) 1 



rn 
N-l 



NT 2 (T 2 - 1 



■ JV-l 



(39) 



where the last line follows from the fact that only m = 
N — 1, N give non-zero contributions to the sum. Com- 
bining all the results, we find that 



N 

i=i 



NT(\T\ 2 - 1)4 (JV-i) 



(40) 



The remaining task is to compute the normalization 
factor Q2i \ a ni\) 2 for T = 1 - 2 1 / N . Recalling the ex- 
pression (|12|l for the matrix elements aki, inserting the 
expression @ for the Jacobi polynomials, and changing 
the indexation such that k and I now range from to 
N — 1, we get 



y\ / / h III 



III 



N 



(T 2 — iy n T l [ " )T- N + )T- k 
\m J \m J 

(41) 

To find the cofactors Am of a = (akl)o<k,KN—i, we make 
use of the fact that they are given by 



A ki = det(a)(a~ 1 ) 



lk- 



We need Afj-ij, and this can thus be found from solving 
the equation ax = 6, where b k — 5k t N-i, and Ajv-i,i = 
det(a)xi. The equation ax = b gives rise to the following 
N equations in x: 



N 

E 

m=0 



V 



T 



-N 



where we have defined the N + 1 quantities y m as 



y m = (T 2 - iy 



El I + 171 

1=0 



T l xi 



(42) 



(43) 



We will now proceed by inverting the two systems l)42|) 
and H43f) in succession. The former system is of course 
underdetermined, which means that the solution for y m 
will contain one free parameter. This parameter will be 
solved for when solving (|43|l for xi . Let us now choose as 
free parameter 77 with 



N 

E 

771=0 



V 



T 



-N„ 



(44) 



then l|42ll turns into 



AT 

E 

777 = 



Vra = {Sk,N-l - 1])T k 



(45) 
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Since k is at most N — 1, the term for m = N in the 
left-hand side vanishes and we can as well sum from to 
N — 1. This system is now easily inverted using the sum 



JV-l 



fe=0 



= 4 



giving 



JV-l / N 



JV-l 



k=0 



(46) 



for < m < N — 1. Noting that the first term in H46|) is 
actually zero for all to < N—l, and T™" 1 for to = N— 1, 
and simplifying further, we find 



r]{T-l) m , 0<to<jV-2, (47) 
N-i _, m -nm ^ 



Dm = 

2/jv-i = T"- 1 - 77(T-l) m . 
The value of un follows from l|44|) as 

JV-l 



Vn 



t n v-y: 



N 



m=0 
N-l 



Vn 



T 



N 



e(:)p , -i)' 



m=0 

T](2T N — (T — l) N ) — NT 



NT 



JV-l 



N-l 



(49) 



We can now solve l|43[l for xi and 77. To perform the 
inversion, we use the sum 



JV-l 



'N-l 



m=0 \p=0 v ' x 

giving, for < I < N — 1, 

N-l /JV-l / N , , 



I + TO 
11) 



m=0 

N-l 



^ w w ; (t 2 - iy 



(N-l 
p=0 

N-l 



I J \mj I (T 2 - l) r 



k=0 



k+m ( m \rpk 



+ (-1) 



i-JV+1 



-rJV-1 



I J (T 2 - l)^- 1 



-^|(?)(^t)' 



+ (-1) 



i-JV+l 



r 



/V-l 



(50) 



We can solve 77 from (|43|l and (|49() . using J5UJ: 
77 [2T* — (T — 1)*] - jVT^" 1 



i=0 



jV 



1=0 

N-l 



p=0 



T 



I \T + 1 



i=0 



jV 



,Y^-^ t jv-i 



S<-*C + /)G)-<-< 



Using the sum 

JV-l 

jv y v^y 

this simplifies to 
rj[2T N ~(T-1) N ] - NT"" 1 



(51) 



- V (T 2 -1) N 
+ {T 2 - 1) N N 



T 



T+l 

rpN-l 



N 



T 



T+l 



JV' 



(T 2 - l)^ 1 
= -r)(T - 1) N {1 - (-T) N ) + N(T 2 - l)T N -\ 

giving, after solving for 77 and simplifying, 

NT 

77 



(2- (1-T) W )' 



(52) 



Inserting this in l|5U|l gives the final solution for x\. 
However, we actually need An-i,u which is given by 
^4-JV-i,f = det(a)x;. With the value for det(a) calculated 
as (O, this becomes 



Ajv-i,i = 

(T 2 - l) iHi ^ [2 - (1 - T) N ] T 
NT 



JV-l 



(2- (1 -T) N ) 



+(-D 



i-JV+1 



ifN-1 



p=0 

jrJV-l 



T 



I VT+1 



I ) (T 2 - l)^- 1 

JV-l 

N(T 2 - l) iUJ 5= il (-l)-'+ 1 T- , + 1 E 

p=0 



T 



I \T+1 



+ (-l) N - l+1 (T 2 - 1) ' ^ ' [2 - (1 - T) w ] 



N- 1 



T 



JV-i-l 



(53) 



We now specialise to the case of interest, namely where 
T takes the optimal value T = 1 — 2 1 / Ar . Then the second 
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term of l|53(l vanishes, giving 



Y 



N (T 2 - l) ini ^ i2 (-l)- i + 1 T- i + 1 ^ 



p=0 



I \T+1 



(54) 



We have not substituted the optimal value for T in every 
part of the expression, because no simplifications occur 
at this point. 

We finally have to add up |j4jv-i,i| over all We first 
show that in the optimal T, the A^-i.i alternate in sign. 
The factor that could change sign in (— l) N ~ l ~ 1 Ajy^i j is 



N-l 



p=0 



I \T+1 



The generating function of this sequence is 

N—l 



S(z) = g i zl = 



i (T + iy-(T + z) 



1=0 



(T + (l-z) 



whence it follows that si can be rewritten as 

N 

J 



si 



I y n )t^ 

(T + l)N-i 2^ \ j ) 

J — £~r -L 

^ N—l— 



t=0 



T*. 



We consider now the functions /j,jv(T) := J2t=o Ct )^ 
and show that they are non- negative for T > — 1/N, 
which covers the case of the optimal T = 1 — 2 1 / jv . First 
put T = z — 1/N , which gives after some simplifications 



1=0 



N\iA fN-l 



KM* - UN) = $>M 7 J e ( t (-W 



This is non-negative if the coefficients 
St=o { N ^ 1 ) (~ are non-negative. This in turn 

would follow from the statement fj-i.N-i(—l/(N — I)) > 
0. Now the terms in the sum 



constitute an alternating sequence, with the terms for 
even t positive and the terms for odd t negative. Non- 
negativity of the sum then follows from the fact that the 
terms decrease in absolute value: indeed, the absolute 
value of the i-term divided by the t — 1 term is (N — t + 
1)/Nt, which does not exceed 1 for t > 0. This finally 
shows non-negativity of the sequence s; . 



As a consequence, we can now just add up all 
(— 1) ~' _1 Ajy_i,i and take the absolute value after- 
wards. This gives 



N-l 



EC-l)"-'- 1 ^-!, 



1=0 



N-l N-l 



N(T 2 - 1)^(-1)-T £ £ Mr' ' 



1=0 p=0 
N-l 



p=0 



N(T 2 - V^^i-l^T J2 (1 + l/TY f j^j 



JV(JV-l) -r 

N 2 T{T 2 -\)^-^{-l) N . 



The final result is then 

2 



( E \ A »™\) = ^ 4 r 2 (T 2 - l)^- 1 ) . (55) 

\m=l / 

Combining the results from Eqs. (|40|l and l|55|) and in- 
serting into the expression H±3|l for the maximal success 
probability we obtain 



1 

7V2 ' 



(56) 



just as we set out to show. This means that there exists 
a beam splitter network with an iV-dimensional ancilla 
state containing photon numbers n; = . . . N — 1 such 
that a generalized nonlinear sign shift can be performed 
on an N + 1-dimcnsional signal state with a probability 
of 1/N 2 . 



IV. CONCLUSIONS AND DISCUSSION 

In this article we have shown how the abstract view 
on linear-optical networks can be used to derive scaling 
laws for success probabilities. Thus far, we have lim- 
ited ourselves to single-shot gates without conditional 
feed-forward dynamics which could in principle be in- 
corporated by concatenating several of those abstract 
networks. We found that the maximal probability of 
success of conditionally generating a (single-shot) gen- 
eralized nonlinear sign shift gate on N + 1-dimcnsional 
signal states using iV-dimensional ancillas with the low- 
est possible photon numbers n; = . . . N — 1 scales as 
1/N 2 . To our knowledge, this is the first time such 
scaling laws have been found. It also hints toward scal- 
ing of success probabilities of certain classes of JV-qubit 
gates. This is due to the fact that multi-qubit quan- 
tum gates acting on tensor-product states with constant 
photon number can be decomposed into a multi-mode 
Mach-Zehnder interferometer-type setup where single- 
mode conditional gates are inserted into the interferome- 
ter's paths 0| . Note, however, that this is not necessarily 
the optimal way to implement such gates. For example, 
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a controlled-fj 2 gate would work in 1/16 of all cases only, 
whereas a more general network has been found in [lq | 
that works with a slightly higher probability of success 
of 2/27. In order to find proper upper bounds on such 
networks, the existing abstract network would have to be 
modified by replacing the single 'active' beam splitter by 
a U(2M)-network, M being the number of modes to be 
acted upon. 

We have shown that computing the success probabil- 
ity within the framework of the abstract network reduces 
to the calculation of various Vandermonde-type determi- 
nants. To do so, we defined a class of polynomials re- 
lated to the Jacobi polynomials. These ^-polynomials 
obey a three-term recursion relation which is given in 
Appendix^ and the elementary symmetric polynomials 
associated with them are derived in Appendix FBI 

We have restricted our attention to the experimentally 
most accessible case in which the ancilla state contains 
the lowest possible photon numbers n; = . . .N— 1. This 
choice represents the only analytically solvable case thus 
far but, on the other hand, is motivated by the fact that 
low photon numbers also means low decoherence which 
is desirable for possible applications in quantum infor- 
mation processing. However, the theory presented above 
is in principle valid for any admissible choice of ancilla 
states. It should be added that the proof technique for 
deriving upper bounds by considering dual convex opti- 
mization problems 01 can similarly applied to the situ- 
ation considered in this article. 
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APPENDIX A: SOME PROPERTIES OF THE 
5-POLYNOMIALS 

In this appendix, we briefly summarize some properties 
of the 5-polynomials. We recall again their definition 
Pity: 

S { f ) {l)=Pf' l - j) {2x 2 -l), 
in terms of the Jacobi polynomials 

P^\x) 

f n + a\ ( n + f3 



2-E 



m I \n — m 



From this definition one can derive a three-term recursion 
relation as 



kS ( k x) (n) = [(x 2 -l)(n + k) + 2k-l]Sl x ! 1 (n) 

-(k-l)x 2 si%(n). (Al) 

Then, according to Favard's Theorem which states 
that if P = (P n ) n >o is a polynomial sequence which sat- 
isfies 

1. P n+1 {x) = {A n x+B n )P n {x)-C n P n - 1 {x) ] P {x) = 
1, P-i{x) := and 

2. A n A n ^C n > 0, 

then there exists a positive measure fi on the real line 
such that P is an orthogonal sequence with orthogonal- 
ity measure /i, the polynomials S^ln) do form an or- 
thogonal sequence for all x ^ 0, ±1. This can be seen by 
noting that k 2 A n A n -iC n — {x 2 — l)x 2 in our case. As 
for the excluded values of the parameters x, we have that 

S { k ±1} (n) = 1 and s[, 0) (n) = it) which both do certainly 
not form orthogonal sequences. Orthogonality is impor- 
tant if one needs to find an inversion relation between 
the powers of n and the polynomials (n) . 



APPENDIX B: SYMMETRIC POLYNOMIALS 
OVER THE S-POLYNOMIAL BASIS 

We now present a number of technical results about 
the expansion of binomial coefficients and related poly- 
nomials in terms of the ^-polynomials. 

For any positive integer N, for any I, and for an integer 
p satisfying p < N, 



N 



with coefficients 



8N - j (x;p;N)= M, 2 ^)(l-x 2 ) w - 



(Bl) 



(B2) 



that are polynomials of degree N — j in x 2 . Specifically, 
for p = N these formulas simplify to 



N 



N 



(x 2 - i) N ( I) = E(-i) w ^' ( N )x 2{N - j) Sf\l). 

3=0 ^ 3 ' 



(B3) 

To prove i|Bl|) and i|B2|) . we work out the right-hand 
side of HBlf) with l|B2|l and plugged in: 



N , \ j / -\ // i 

j=0 



r=0 



T \ r 



(x 2 -iy 



Since the binomial coefficients take on the value for 
j < 0, j > p, r < and r > j, and since we have 
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the condition p < N, we can drop the bounds on the 
summation signs and subsequently move them up front, 
giving: 



Now we notice (p (^,) = ( p ) , so that we can calcu- 
late the sum over j easily, using 

i(_x 2 )^ = (1 -x 2 ) 



J-r 



(P M C_~2\p-j - n _ ^Np-r 



It remains to calculate 



E Q (' t r ) (1 _ ^""V - i) w - p+r , 

which is nothing but 



It now remains to show that the sum in this expression 
reduces to (^). To do that, we note that (^) is a poly- 
nomial of degree p in I, with zeroes 0, 1, ... ,p — 1 and 
leading term l p /pl. We will confirm that the same holds 
for the sum J2 r CO CV) (-!) p ~ r - Consider first the spe- 
cific values I — 0,1, ... ,p— 1. The factor ( 1] r ) is a degree 
I polynomial in r. It can therefore be written as a linear 
combination of the degree j polynomials (p . Thus, with 
cj independent of r, 



p 

L 

r=0 



ST lP\f l + r \( -np-r 



(-iy 



r I \ r 



I r 



3=0 
I 



r J \3 



(~iy- r 



^(of!)E 



J / r=0 V J 



P- 3 



{~l) p - r . 



The sum over r in the latter expression is clearly for 
j < p, hence the whole expression is indeed zero for I = 
0,1, . . . ,p — 1. Now, concerning the sum's leading term 
in I, we note that it can only come from the term with 
summation index r = p, ( , whose leading term in I is 
l p /pl, as required. This completes the proof of l|B2|) . 



The binomial coefficient ( J is a polynomial in I: 



aw-i). ..(i-p+i) = ^na-fe) 



PJ p 



k=0 



We now look at related polynomials where one or two of 
the factors (I — k) are missing. 



Consider first the polynomial 

i p i 



l-q pi 



fc=0 
kjtg 



where q is an integer in the range < q < p. This 
polynomial can be expanded in terms of the binomial 



coc: 



efficients Q , with r < p: 



p-i 



1 'i ) v//v 



The proof of this expansion goes as follows. Consider 
the numerator: 



p-i 



E(-d 



-l-r f r \ ( l 



r=0 



qj \r 



This can again be rewritten using 
p-i 

i-iy- 1 -^ 



±\-iy-^-°( l - q 



s=0 

l-q-1 
p-q-1 



Wrapping up everything yields the expression 
(g) (p-^-i) 



which is very easily seen to be identical to (^) /(I — q). 

Using this expansion and l|B3|) , one can easily show the 
following. For any positive integer N , for any I, and for 
an integer q satisfying < q < N, 

( I \ N 

{x ^SimL =J2(-l) N - j 4l .(z;JV)sf (Z), (B4) 
l q 



3=0 



with coefficients 



,(?) 



>(x;N) 



1 " 



.N—r 



N 



\ r W r \ x 2{r ~ i) {l - x 2 ) N ~ r 



(B5) 
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APPENDIX C: SIMPLE VANDERMONDIANS 
WITH GAPS IN THEIR ARGUMENTS 

We have seen that we need to calculate simple Van- 
dermondians of degree N — 1 with integer arguments 
[0...N - l]\{k} [see Eq. ©]. Because of Eq. 
it is enough to concentrate on Vandermondians over 
the power basis. First, note that the Vandermondian 
Vn(xi, . . . , xn) can be decomposed as 

JV-l 

V N (xi, . . .,x N ) = JJ (x N - x t )V N -i(xi, . . . ,x N -i) . 



s-l N-s 

YlvRp =(*-mN-s)\. (C2) 

9=1 P=l 



N 

Noting that V N = ]J (n-1)! = Vjv_i(JV-1)!, we obtain 

71=1 



Hence, we can write 



V, 



(8-1) 



iV-l 



The product in the denominator can be computed as 



(CI) 



N 



N 



n \ s -™\= 

m=l, 771.7-s m— 1 



= [| (s - to) (to - s) 



m— s+1 



(C3) 



With fc = s — 1 and the additional factor (\T\ 2 

1) (iV-l)(iV-2)/2 ) thig ig just Eq ^ 
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